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ABSTRACT 



/K\ms. We study the variable star content of the globular cluster M53 to compute the physical parameters of the constituting stars and 
the distance of the cluster. 

Methods. Covering two adjacent seasons in 2007 and 2008, new photometric data are gathered for 3048 objects in the field of M53. 
By using the OIS (Optimal Image Subtraction) method and subsequently TFA (Trend Filtering Algorithm), we search for variables 
in the full sample by using Discrete Fourier Transformation and Box-fitting Least Squares method. We select variables based on the 
statistics related to these methods combined with visual inspections. 

Results. We identified 12 new variables (2 RR Lyrae stars, 7 short periodic stars - 3 of them are SX Phe stars - and 3 long-period 
variables). No eclipsing binaries were found in the present sample. Except for the 3 (hitherto unknown) Blazhko RR Lyrae (two RRab 
and an RRc) stars, no multiperiodic variables were found. We showed that after proper period shift, the PLC (period-luminosity- 
color) relation for the first overtone RR Lyrae sample tightly follows the one spanned by the fundamental stars. Furthermore, the 
slope is in agreement with the one derived from other clusters. Based on the earlier Baade-Wesselink calibration of the PLC relations, 
the derived reddening-free distance modulus of M53 is 16.31 ± 0.04 mag, corresponding to a distance modulus of 18.5 mag for the 
Large Magellanic Cloud. From the Fourier parameters of the RRab stars we obtained an average iron abundance of - 1 .58 ± 0.03 (error 
of the mean). This is ~ 0.5 dex higher than the overall abundance of the giants as given in the literature and derived in this paper from 
the three-color photometry of giants. We suspect that the source of this discrepancy (observable also in other, low-metallicity clusters) 
is the want of sufficient number of low-metallicity objects in the calibrating sample of the Fourier method. 
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1. Introduction 

Among the oldest and most metal-poor Galactic globular clus- 
ters, the outer halo cluster M53 (NGC5024, a = 13*'12'"55!3, 
6 = -h18° 10'9") is the second most ab undant in variable stars af- 
terM15 jClement et al.ll200U see also lMarfn-Franch et al.ll2009l 
for cuiTent age deter mination s ). The search for its RR Lyrae 
content dates back to ' ShaplevI (Il920h . and yielded a total of 60 
RR Lyrae stars until now. A thorou gh review of discoveries is 
given in the paper of lKopackH (|2000|) . the only CCD time-series 
photometric study of M53 RR Lyrae stars prior to this work. 
This cluster is also very abundant in blue straggler stars (BSSs). 
Amon g the al most 200 BSSs, 8 were pro ven to be SX Phe vari- 
ables jjeon et a l. 2003). We also note that lBeccari et all (l2008al) 
showed that the BSSs follow a bimodal distribution, implying a 
binary rate of 10%, but no eclipsing binaries have been found 
so far As for t he most extensive color-magnitude study of M53, 
we refer to Rev et al.l (11998 ). The cluster is located at a very 
high Galactic latitude of b = 79°76. Therefore, its field contam- 
ination and interstell ar reddening is expected to be negligible. 
This is confirmed by Zinn' ('1985'), reporting E(B - V) = 0.0, 
and by Schlegel et al. (1998), who give a reddening value of 
E{B - y) = 0.02. Despite its high variable content and its fa- 
vorable position, no comprehensive wide-field CCD time-series 
photometric study has been publishe d for this cluster. Among 
the two previous works of similar type, Kopacki (2000) observed 
only the central part of the cluster, while Jeon et al. (2003) lim- 
ited themselves to the investigation of SX Phe stars. The present 
study, based on two-color photometry covering most of the clus- 



ter, intends to improve the time and spatial coverage of this clus- 
ter. 

To exploit the entirety of our photometric data-base and 
detect as many variables as possible, we employ the widely 
used method of O ptimal Image Subtractio n (OIS ) devel oped by 
lAlard & LuptonI ( Il998.) and improved by AlardI (l2000h . to ob- 
tain optimum relative light curves. Then, in a post-processing 
phase, we filter out temporal trends by using the Trend Filtering 
Algorithm (TFA) of Kovacs et al. (2005). This latter method 
has been extensively used during the past several years by 
the HATNet ProjecQ in search for transiting extrasolar plan- 
ets and also employed successfully in gene ral variable search 
(iKovacs & Bakosl2008t[Szul agvi et al.'2009). The main purpose 
of this paper is to search for variables in the full dataset and not 
limit the investigation to the well-known regions of variability 
(i.e., the RR Lyrae strip, the blue straggler and the upper red 
giant branch regions of the Hertzspring-Russell diagram). With 
the aid of the reduction and post-processing methods mentioned 
above, we can increase photometric precision and attempt to 
reach the detection limit set by the photon statistics. 

In Sect. |2] we give details on the observations and data pro- 
cessing, and show examples of the detections resulting from the 
methods applied. The next three sections (Sects. |3}|5]l are de- 
voted to the more detailed analysis of the specific types of vari- 
ables, most importantly to that of the RR Lyrae stars. Our at- 
tention will be focused mostly on the period-luminosity-color 



' Hungarian-made Automated Telescope Network, see 
http : / /cf a-www . harvard . edu/~gbakos/HAT 
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relation and the metallicity issue. A brief discussion and conclu- 
sions are given in Sect. |6] 

2. Observations, data reduction, method of analysis 

2.1. Observations and data reduction 

We performed time-series photometric observations of M53 
through V and / filters of the Johnson and Kron-Cousins sys- 
tem, respectively, with the 60/90/180cm Schmidt Camera of 
the Konkoly Observatory located at Piszkesteto. The telescope 
is equipped with a Photometries AT 200 type CCD of 9-fi pix- 
els in a 1536 x 1024 array. The point spread function (PSF) is 
slightly undersampled with a scale of 1.026 arcsec/pixel. The 
observations were carried out between March 25, 2007 and May 
28, 2008 in two seasons, on a total number of 26 nights. The 
observations have a full time span of 430 days. Exposure times 
varied between 550 and 750 seconds in V and 350 and 540 sec- 
onds in / depending on atmospheric conditions. A total number 
of 340 and 230 frames have been collected in the V and / bands, 
respectively. While for the time series analysis we use only the 
V-band data, for standard magnitude transformations, variable 
star classification, and for deriving cluster distance and metallic- 
ity, we also utilize the /-band data. 

The processing of the images (bias, dark, and flat-field cor- 
rections) and the alignment of all frames into a common pixel 
reference system were performed by standard irai0 packages. 
The joint field of view of all observations was somewhat trun- 
cated compared to the full detector area due to technical rea- 
sons. Only this 13'5 x 14.'5 common area (still covering a fairly 
large part of t he cluster by spanning to ~ 0.7 tidal radius, see 
the catalog of iHarrisll 19961 revised in February 2003), centered 
approximately on the cluster core, was subjected for further re- 
duction. 

We calibrated the astrometric transformations between pix- 
els and celestial coordinates by using 30 bright stars from 
the second editio n of the Guide Star Catalog (version 2.3.2, 
iLasker et al.ll2008l) uniformly distributed around the cluster cen- 
ter and located sufficiently well outside the cluster core. The 
standard deviations in the residuals of the coordinate mapping 
were 0'.'148 and 0'.'163 in right ascension and declination, re- 
spectively. 

In order to deal with heavy crowding and to probe the regions 
close t o the center of the cluster, we applied the OIS method of 
lAlardI (120001) . For this method one has to prepare an accurate ref- 
erence frame (RF) in each band which will be subtracted from 
all other frames after the appropriate scaling of the PSF of the 
stellar objects. The success of the method is highly sensitive to 
the quality of the RFs, therefore we carefully selected the best 
images with the lowest background, best seeing, and most reg- 
ular PSFs and stacked them into one single image in order to 
increase their signal-to-noise ratio. Before the subtraction from 
each frame, the RF was convolved with a kernel function, which 
was determined by using the representative PSFs of 225 stars 
uniformly distributed on the image. The kernel functions were 
allowed to have second order spatial variations to account for a 
varying PSF across the image. The convolved RF, being the best 
approximation of the processed frame in a least-squares sense, 
was subtracted from the processed frame. The method yielded a 



^ IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation 



set of residual images consisting of, in principle (but not in prac- 
tice), the true light variations plus noise. For all procedures of 
ir nage convolut ion and subtraction we used the isis 2.2 package 
of lAlardl(l2000h . 

We applie d the daophot PSF-fitting software package 
dStetsonl 19871) on the RFs to identify stellar objects and to obtain 
their instrumental magnitudes. We found 3048 stars above the 
lOcr level. Linear transformation equations of these instrumen- 
tal magnitudes into standard ones wer e esta blished by using 110 
photometric standard stars of lStetsoiil (I200Q) around the cluster 
center 

We measured differential flux variations on the residual 
frames by aperture photometry at all stellar positions obtained 
from the RFs, using a custom IRAF routine. The sizes of the 
apertures were variable in order to conform to the FWHMs of 
the PSFs on the original frames, changing in time due to see- 
ing variations. To diagnose if some faint variables might have 
been remained undetected on the RFs due to crowding but had 
light variations emerged on the residuals, we also employed the 
detection algorithms of isis 2.2 in two ways. First, all sources 
of significant variability were checked for the case if light vari- 
ation comes from a yet undetected faint star falling within the 
same aperture. In these cases we performed its photometry us- 
ing a new, accurately recentered aperture. Secondly, we checked 
if there were any remaining faint variables that had neither been 
detected on the RFs nor fell within any aperture (we note that 
we did not find any additional variables by this latter test). All 
light curves obtained in this way, expressed in linear flux units, 
were then further processed with the aid of the Trend Filtering 
Algorithm (TFA, see Sect. l2.2l for further details). 

As the last step of the reduction, to convert differential fluxes 
of an object into magnitudes, one has to know the flux /rf of the 
source on the RF in each band which has been subtracted from 
each processed frame (see e.g.. lWozniakll2000l) . This additional 
information sets the zero point of the logarithmic transformation 

Am; = -2.5 log(/RF -1- Afi) 

between the Am,- diff'erential magnitude and A/j differential flux 
of the object (derived by the OIS method) on the i-th frame. The 
/rf value has been measured by PSF photometry for each ob- 
ject as described above. Clearly, a flux-to-magnitude transforma- 
tion with an imprecise /rf value yields an incorrect amplitude of 
the light variation and can significantly alter the magnitude av- 
erages as well, but its actual effect on these errors is specific 
to each object and depends on the brightness and the shape of 
the light variation. Therefore, the flux curves of variable stars 
were transformed into magnitudes only if their (already stan- 
dardized) photometry on the RFs were reliable enough, which 
has been decided based on the formal errors of the PSF photom- 
etry. Generally, objects were transformed into magnitudes if the 
formal error of their (V - /)rf was below 0.03 mag (see Sect. 13.21 
for details on this cut in the case of RR Lyrae stars 

2.2. Post-processing witli TFA 

The time-series data-base obtained above is subjected to a 
subsequent post-processing phase, when we filter out various 
trends/systematics due to, e.g., imperfect reduction, lack of cor- 
rection to position- and time-dependent extinction, anomalies in 

^ Photometric time-series of the variables detected in this work and 
the Fourier decompositions of the RR Lyrae stars with magnitude- 
transformed light curves in V color are available electronically online 
at CDS. 
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Fig. 1. Histogram of the peak frequencies of the Fourier (DFT) 
spectra for the 3048 stars analyzed. Upper panel shows the fre- 
quency distribution of data obtained by aperture photometry on 
the residual frames of the OIS method. Lower panel shows the 
effect of TEA on these data. 



the convolution procedure prior to image subt r action , etc. The 
method is described in detail by iKovacs et al.l (l2005l) and also 
summ arized recently by iBakos et al.l ( 2009 ) and Szulagvi et al.l 
(l2009l) . Here we only briefly note that the method is based on 
the idea of correcting elements of the systematic variation in the 
target time series by using the light curves of many other ob- 
jects, available in the CCD frame. In the course of signal search, 
we have no information on the temporal content of the time- 
series. Therefore, we assume that the target is trend- and noise- 
dominated. If it is the case, then the method suppresses the trend 
and gives rise to the signal to appear in some time-frequency 
transformation of the detrended signal. Here, on the expense of 
the suppression of the systematics, we also deform the true sig- 
nal at some level. This effect is cured in the second step of filter- 
ing, when we reconstruct the signal by using a full time-series 
model that includes both the systematics and the signal with the 
period obtained in the first step. An extension of the me thod to 
multiperiodic signals is given in lKovacs & Bakosl(l2008l) . In this 
method (adopted also in here) one fits the TFA template and the 
signal (represented by a Fourier series) simultaneously, thereby 
avoiding iteration by incomplete time-series model representa- 
tions. 

Special attention is to be paid here to the effects related to 
the number of data points. In ideal case, the necessary num- 
ber of TFA templates to handle all systematics is considerably 
lower than the number of data points constituting the time-series. 
Unfortunately, we do not have a good method to select a "rep- 
resentative" small set of templates, therefore we use a large set 
that includes both "useful" and "useless" (e.g., pure noise) tem- 
plate time-series. Attempts have been made to determine an op- 
timum number of templates but it seems that it is hard to set - 
tle at a value lower than a few hundred (Szulagvi et alj|2009l) . 
Since in the present case we have only 340 data points per time- 
series, we have to limit the number of templates considerably 
(i.e., about half of the number of the data points) to avoid over- 
fitting and sudden increase of the false alarm probability (FAP). 
This latter effect is thoroughly tested in each interesting case by 
using multiple tem plate runs, wh ich are capable of reducing FAP 
by several factors (IKovacs & Bak os 2007). To avoid significant 
overfitting, we settled on the lowest TFA template numbers that 



yielded the most significant detection together with the lowest 
unbiased scatter around the signal. This optimum template num- 
ber varied between 60 and 140, the latter being an upper limit 
set by us due to the low number of available data points. 

In general, we selected the TFA template sets from the 
~ 1000 brightest stars in the sample. The faintest stars in the 
template sets did not exceed ~ 18 mag in V band. The templates 
cover the frame in a quasi-uniform manner (see I Kovacs et al.l 
i200l for more details on the distribution of templates). 

The signi ficance of the Fourier or BLS components (see 
IKovacs etal]r2002) was deduced by checking the SNR of the 
frequency spectra (the ratio of the amplitude of the highest 
peak to the standard deviation of the spectrum). Statistical 
tes ts (similar to th e ones performed in our earlier papers, e.g., 
in iNagv & Kovacsl 12006 ) have led to the conclusion that for 
Gaussian signals with SNR > 6.7, the FAP due to uncorre- 
cted random noise is less than 1%. Considering that the pro- 
cessed signals always have some colored noise (either physical 
or instrumental) we selected variables well above this SNR value 
(usually with SNR exceeding ~ 8). 

In demonstrating the ability of signal detection, first we com- 
pare the distributions of the peak frequency components ob- 
tained by OIS and by the subsequent application of TFA on the 
OIS time series. Figure 1 shows that already the OIS data are rea- 
sonably free of periodic systematics. Actually, when compared 
with earlier similar diagrams derived on data of apert ure pho- 
tometry (e.g., Kovacs et al. 2005; Szulagyi et al.ll2009|) we see 
that the current data are much less dominated by the custom- 
ary 1 d"' systematics. We see in the bottom panel that TFA has 
successfully filtered out the daily trends. In the final dataset the 
frequency distribution becomes nearly flat, with a slight surplus 
due to the RR Lyrae stars and the long-term irregular light vari- 
ations of red giants. 

In checking the effect of TFA in more detail, we plotted the 
DFT (Discrete Fourier Transformation) spectra and the folded 
light curves for the OIS and OIS+TFA data in the case of an 
RRc (V71, Fig.© and an SX Phe star (V78, Fig. [3). The eff^ect 
of the signal reconstruction is large in both cases. The true vari- 
ability of V71 would have been significantly more troublesome 
to decipher from the forest of high peaks in the original (OIS) 
data than in the clean, single-component spectrum of the TFA- 
filtered time-series. For V78 the situation is better because of the 
high frequency of the pulsation. Nevertheless, even in this case, 
the highest peak is near 1 d ' (and its aliases), therefore, an auto- 
matic search that only looks for the highest peak in the spectra, 
would not find this variable in the original (OIS) data. The im- 
provement in the quality of the folded light curve derived from 
the reconstructed data is also substantial. We should note how- 
ever that the small scatter seen is somewhat biased, because of 
the large number of TFA templates used. The effect of the bias in 
the residual scatter is common to all data regressions. The fit in- 
troduces correlation among the original array elements that leads 
to smaller observed scatter as it were expected if only the deter- 
ministic (model-predicted) part of the signal would have been fit- 
ted. In general, when least-squares fitting a data array of data 
points with A^,, parameters (larger than exact model represen- 
tation required) the residual RMS will be -sJ{N - Np)/N-times 
of the v alue expected from th e noise level in the original data 
(see also lSzulagyi et al]|2009l for further details on this subject). 
For this reason, the unbiased standard deviations of the recon- 
structed light curves are 10% and 30% higher due to the 60 and 
140 templates used for V71 and V78, respectively. The different 
template numbers resulted from the multi-template tests men- 
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Fig. 2. Detection of the RRc vaiiable V71. Top left: folded light 
curve (with the period detected by TFA) of the OIS data; top 
right: DFT of the OIS data (light shade) and that of the TFAd 
data (dark shade); bottom left: folded light curve of the TFA- 
reconstructed data; bottom right: DFT of the TFA-reconstructed 
data. Star identification name and period (in days) are shown in 
the upper right corner Please check note in text on the scatter of 
the TFA-filtered light curve. 
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Fig. 3. Detection of the SX Phe star V78. Notation is the same 
as in Fig.|2l 



tioned earlier in this section, in the description of the method of 
analysis. 

In a final example, in Fig.|4]we show the signal reconstruc- 
tion capability of TFA. We used an 8th order F ourier series with 
144 T FA templates in the reconstruction (see iKovacs & BakosI 
12008 ^). Inspection of the DFT spectrum of the OIS light curve 
of this star shows that it is not influenced by periodic systemat- 
ics of type integer d Therefore, the improvement introduced 
by TFA is attributed to filtering out some transient signals. We 
found transient syste matics also in o ther v ariables, similarly to 
the ones reported by Szulagvi et al.l (l2009t) . The unbiased stan- 
dard deviation of the residuals around the best-fitting Fourier- 
sum dropped down by 14%, which is a good improvement over 
the original data. 



V45 



0.6549513 




0.2 



0.4 0.6 

phase 



0.8 



Fig. 4. Signal reconstruction in the case of the RRab star V45. 
Open circles denote the OIS fluxes while filled ones show the 
TFA-processed values. Please check note in text on the scatter 
of the filtered light curve. 



3. RR Lyrae stars 

Our prime interest in this paper is the study of the RR Lyrae stars, 
since they are in large number in the cluster and they are the ones 
that have fairly well established theoretical and observational un- 
derstanding (e.g., variable and mode identification, computation 
of the physical parameters from observed ones, etc). We focus 
on the distance and on the metallicity (i.e., [Fe/H]) as the two 
most important parameters derivable from the available empiri- 
cal relations. 



3. 1 . General description 

We identified altogether 54 RR Lyrae stars in the cluster Except 
for the two RRc stars V71 and V72, all others were known pre- 
viously (see the C atalog of Variable Sta rs in Globular Clusters, 
hereafter CVSGC, IClement et all 12001') . We note that the vari- 
ability status of V34, which had be en questioned previously by 
Ivan den Hoven van GenderenI (1 19471) was clearly confirmed by 
our data. Due to the application of OIS and TFA, many of the 
variables ended up with light curves suitable for further analy- 
sis. We give the full list of variables together with their basic 
observed parameters in Table [T] Please note the very high SNR 
attached to each detection even for the variables sitting near the 
center of the cluster 

The accuracy of the periods of the RR Lyrae stars has been 
estimated in the following way. First we computed the Fourier 
decompositions from the TFA-reconstructed light curves. From 
these decompositions synthetic light curves were generated and 
Gaussian noise was added with standard deviations obtained 
from the unbiased estimates of the residual scatter between the 
synthetic light curves and the TFA reconstructions. The gen- 
erated light curves were subjected to the standard Fourier fre- 
quency analysis as if they were real observed Ught curves. The 
periods obtained from each realization were stored and after 
completing 100 independent simulations for each object, the 
standard deviations of those period values were computed. The 
errors shown in Table[T]are these standard deviations. 

The periods obtained in this study are generally in good 
agreement with previously determined values, except for the 
three RRc stars V44, V47, and V55, for which previous val- 
ues are various aliases of the true periods. In the case of other 
first overtone stars the agreement between our and earlier pub- 
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Table 2. Modulation properties of the Blazhko-stars 



ID 


Type 




Pm [d] 




Vll 


RRab 


0.028 


35 


0.21 


V16 


RRc 


0.017 


58 


0.41 


V57 


RRab 


0.020 


52 


0.25 



Note: j?m 1 is the ratio of the amplitude of the largest modulation peak in the 
DFT spectrum and the amplitude of the /o pulsation component (see Table[T), 
in V flux. 



lished values are better than 0.0004 d ' with a typical difference 
of about 0.00003 d As for RRab stars, deviations are gener- 
ally below 0.00002 d"\ with a few higher values not exceeding 
0.00012 d-i. 

Figures|5]and|6]show the TFA-reconstructed and folded light 
curves. Standardizing the photometry for the highly crowded 30 
objects shown in Fig. |6] was not reliable enough, therefore, we 
left the light curves in relative flux units. Unfortunately, these 
objects can only be partially utilized here due to the lack of em- 
pirical relations derived on fluxes and the absence of existing 
photometry with a higher resolution CCD (exceptions are V7, 
V8, V9, V33, and V37 tha t have good-quality light curves in V 
magnitude from .Kopackil 12000) . We note that among the pre- 
viously known RR Lyrae stars missing from our sample, V52 
and V53 were merged with each other, V61 was merged with 
the long period variable V49, and the rest (V12, V13, V21, V30, 
and V48) were outside the field of view. 

3.1.1. Blazhko-stars 

We investigated possible multiperiodicity by performing succes- 
sive prewhitening on all 3048 objects, including RR Lyrae stars. 
Surprisingly, from the whole sample we found only the three 
Blazhko stars given in Tables [1] and |2] in spite of the relatively 
large scatter visible in several stars (see Fig.|6]l. We note that the 
detection of the Blazhko behavior in the three stars is firm. The 
modulation is clearly observable as side-lobe frequencies with 
uniform /,„ frequency separations after prewhitening the main 
frequency and its harmonics. For Vll, significant side lobe fre- 
quencies could be identified for up to the 7th harmonic of /o at 
n/o + fm, except for n - 2, where we found 2/o - /m. In the 
case of the RRc star V16 we found /o + and 2/o + /m. The 
phase modulations in both stars are very little, unlike to that of 
the third Blazhko-star, V57, whose modulation is very large both 
in phase and in amplitude. We detected /o +fm and 2/o -H/m in the 
case of V57. The light curves of the Blazhko stars folded with 
the pulsation periods are shown in Fig.|6l The light curves were 
reconstructed by simultaneously fitting the template and the full 
modulated signal (see Sect. 12.21 ). 

3.2. PLC relation 

The tightest empirical relations existing for RR Lyrae stars are 
the period-infrared K magnitude (PL^:) (' Longmore et al.l[l990l) 
and the various period-luminosity-color (P LC) relations (see 
Kovacs & Walker 2001) fherea fter KWOl]; ICortes & CatelanI 
2008- .Caceres & CatelanI pOOSh. In addition to the theoretical 
interest in these relations dBono et a n i200lH D i Criscienzo et alj 
I2OO4 I). they have also important practical utilization by giving a 
simple and accurate way to determine distances. Furthermore, it 
is also an interesting question whether RRc stars (after a proper 
period shift) fit the empirical PLC relation spanned by the RRab 
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Fig. 7. PLC relation for 24 RR Lyrae stars in M53. RRc periods 
were shifted by adding 0.128 to their log Pi. Continuous line 
shows the zero-point fitted regression line with the slope given 
in KWOl. 



stars. Since M53 contains a considerable number of RRc stars, 
we can a ddress this qu estion with reasonable confidence here 
(see also ICassisi et~aL I [2004 for concerning the same problem 
in the case of the variables in M3). To decrease the number 
of parameters fitted, we use the slope of the PLC relation (i.e., 
log P — Wesenheit index relatiorQ) derived by KWO 1 from glob- 
ular cluster RRab variables, for the V and / magnitudes. (We 
note that this slope is in agreement with the one derived from the 
large - but considerab ly noisier - sample of LMC RRab stars by 
ISoszvnski et"ani2003l ) 

Figure |7] shows the derived joint PLC relation. The periods 
of the RRc stars have been increased by -log(Pi/Po) - 0.128 
to take proper account for the overall period difference be- 
tween the RRab and RRc stars. The above shift corresponds to 
Pi/Po = 0.745, which is the average expected period ratio for 
the first overto ne and fundamen tal mode periods for RR Lyrae 
stars (see, e.g.. ICox et al.|[T983l) . From the 54 RR Lyrae stars 
with available V, I photometry we omitted 30 stars due to their 
outlier status in the PLC relation (the stars without magnitude 
averages in Table [TJ. We investigated the question if the omis- 
sion of all these stars can be justified on the behavior of some 
other parameters derived during the reduction of the data. The 
first parameter to be utilized for this purpose could be the scatter 
of the light curve. But while crowding introduces large errors in 
the zero points of the flux-to-magnitude conversion, this may not 
necessarily lead to the increase of the scatter of the light curve 
due to the effectiveness of the OIS method (compare variable 
positions in Table[T]and light curve scatters in Fig|6]l. Therefore, 
we resorted the inspection of another quantity that might corre- 
late with the zero point errors of the light curves. This quantity 
was the formal errors of the RF magnitudes as they came out 
from the standard method of PSF photometry. For a consider- 



By construction, the Wesenheit index is redd ening-free in the cas e 
of "standard" interstellar reddening law (e.g., ICardelli et alj Il989h . 
Interestingly, the best-fitting PLC relation with free-floating color co- 
efficients yield v alues very close to those predicted by the standa rd ex- 
tinction law (see lKovacs & Jurcsi]dri997l ; and iUdalski et al.lll999l) . 
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Table 1. Properties of the RR Lyrae stars 



ID 


Type 


a 


[hms] 


S [dms] 


dl"] 


Pld] 


SNR 


V 


V-I 


VI 


RRab 


13 


12:56.3 


18:07:13.8 


176.8 


0.6098298(9) 


14.3 


16.970 


0.579 


V2 


RRc 


13 


12:50.3 


18:07:00.8 


202.0 


0.386122(1) 


22.2 


16.879 


0.402 


V3 


RRab 


13 


12:51.4 


18:07:45.5 


154.6 


0.630605(1) 


17.8 


16.848 


0.447 


V4 


RRc 


13 


12:43.9 


18:07:26.4 


230.3 


0.385545(1) 


26.0 


16.800 


0.383 


V5 


RRab 


13 


12:39.1 


18:05:42.6 


353.3 


0.639426(1) 


23.2 


16.888 


0.526 


V6 


RRab 


13 


13:03.9 


18:10:19.8 


123.7 


0.664020(1) 


14.4 


16.819 


0.522 


V7 


RRab 


13 


13:00.9 


18:11:29.7 


113.3 


0.5448584(6) 


13.5 






V8 


RRab 


13 


13:00.4 


18:11:05.1 


92.3 


0.615528(1) 


16.1 






V9 


RRab 


13 


13:00.1 


18:09:25.0 


81.9 


0.6003690(7) 


15.2 






VIO 


RRab 


13 


12:45.7 


18:10:55.6 


143.3 


0.6082612(7) 


18.8 


16.837 


0.463 


Vll 


RRabB 


13 


12:45.4 


18:09:01.9 


156.3 


0.629940(5) 


12.2 






V14 


RRab 


13 


13:20.5 


18:06:42.8 


415.6 


0.5454625(7) 


14.0 


16.880 


0.452 


V15 


RRc 


13 


13:12.4 


18:13:55.0 


332.5 


0.3086646(9) 


21.6 


16.894 


0.342 


V16 


RRcB 


13 


12:46.2 


18:06:39.2 


247.3 


0.3031686(7) 


18.7 






V17 


RRc 


13 


12:40.4 


18:11:54.1 


236.7 


0.381282(1) 


27.5 


16.848 


0.422 


V18 


RRc 


13 


12:48.7 


18:10:12.9 


93.9 


0.336054(1) 


22.2 






V19 


RRc 


13 


13:07.0 


18:09:26.1 


173.2 


0.391377(1) 


28.1 


16.880 


0.437 


V20 


RRc 


13 


13:09.0 


18:04:16.2 


404.7 


0.384337(1) 


22.5 


16.875 


0.407 


V23 


RRc 


13 


13:02.3 


18:08:35.9 


137.9 


0.365804(1) 


25.2 


16.825 


0.391 


V24 


RRab 


13 


12:47.2 


18:09:33.0 


120.4 


0.763198(2) 


20.3 






V25 


RRab 


13 


13:04.4 


18:10:37.2 


133.5 


0.705162(1) 


16.7 


16.779 


0.555 


V26 


RRc 


13 


12:35.7 


18:05:20.5 


401.7 


0.391106(1) 


24.1 


16.845 


0.395 


V27 


RRab 


13 


12:41.4 


18:07:23.9 


257.9 


0.671071(1) 


14.7 


16.856 


0.550 


V28 


RRab 


13 


12:42.1 


18:16:37.9 


430.9 


0.6327804(7) 


16.5 


16.875 


0.545 


V29 


RRab 


13 


13:04.3 


18:08:46.9 


152.9 


0.823243(4) 


26.6 


16.782 


0.588 


V31 


RRab 


13 


12:59.6 


18:10:04.6 


61.9 


0.705665(1) 


15.2 






V32 


RRc 


13 


12:47.7 


18:08:35.9 


142.6 


0.390623(2) 


20.4 


16.712 


0.416 


V33 


RRab 


13 


12:43.9 


18:10:13.2 


162.5 


0.6245815(8) 


16.9 






V34 


RRc 


13 


12:45.7 


18:06:26.1 


262.0 


0.289611(1) 


21.2 


16.928 


0.314 


V35 


RRc 


13 


13:02.4 


18:12:37.5 


179.8 


0.372666(2) 


22.7 


16.943 


0.442 


V36 


RRc 


13 


13:03.3 


18:15:10.2 


321.8 


0.373242(1) 


21.5 


16.879 


0.414 


V37 


RRab 


13 


12:52.3 


18:11:05.1 


69.6 


0.717615(1) 


21.0 






V38 


RRab 


13 


12:57.1 


18:07:40.5 


151.8 


0.705792(2) 


21.2 


16.749 


0.557 


V40 


RRc 


13 


12:55.9 


18:11:54.7 


105.3 


0.3147939(9) 


19.6 






V41 


RRab 


13 


12:56.8 


18:11:09.3 


63.6 


0.614438(1) 


19.2 






V42 


RRab 


13 


12:50.6 


18:10:19.5 


67.4 


0.713717(2) 


21.8 






V43 


RRab 


13 


12:53.1 


18:10:55.7 


55.5 


0.712017(2) 


24.2 






V44 


RRc 


13 


12:51.5 


18:09:58.2 


54.6 


0.375099(2) 


23.7 






V45 


RRab 


13 


12:55.2 


18:09:27.4 


42.4 


0.654950(2) 


16.2 






V46 


RRab 


13 


12:54.6 


18:10:36.2 


28.1 


0.703655(3) 


22.0 






V47 


RRc 


13 


12:50.4 


18:12:24.6 


151.4 


0.335377(1) 


21.1 


16.791 


0.372 


V51 


RRc 


13 


12:57.8 


18:10:50.7 


54.7 


0.355203(2) 


20.2 






V54 


RRc 


13 


12:54.4 


18:10:31.4 


25.0 


0.315122(3) 


14.2 






V55 


RRc 


13 


12:53.6 


18:10:39.2 


37.7 


0.443386(2) 


26.4 






V56 


RRc 


13 


12:53.7 


18:09:29.1 


46.3 


0.328796(2) 


15.1 






V57 


RRabB 


13 


12:55.5 


18:09:58.0 


12.1 


0.568234(7) 


11.8 






V58 


RRc 


13 


12:55.6 


18:09:30.8 


39.4 


0.354954(3) 


17.2 






V59 


RRc 


13 


12:56.7 


18:09:20.4 


53.4 


0.303941(2) 


19.7 






V60 


RRab 


13 


12:57.0 


18:09:36.0 


41.6 


0.644755(2) 


18.4 






V62 


RRc 


13 


12:54.0 


18:10:30.1 


27.0 


0.359891(4) 


13.8 






V63 


RRc 


13 


12:56.2 


18:10:02.8 


15.5 


0.310476(4) 


14.3 






V64 


RRc 


13 


12:52.6 


18:10:12.5 


38.5 


0.319529(1) 


19.2 






V71* 


RRc 


13 


12:54.5 


18:09:54.3 


18.8 


0.304242(5) 


14.5 






V72* 


RRc 


13 


12:55.8 


18:09:50.6 


20.7 


0.254155(3) 


13.5 







Notes: 

- The radial distance d ["] is measured from the nominal cluster center defined in Stetson's catalog of photometric standards l lStetsonir200(]| : 
a = 13''12'"55'.26; S = +18°10'09'.'8 [2000.0]). 

- Errors corresponding to the last digits of the periods are shown in parentheses (see Sect. l3.1l for details on the computation of the errors). 

- The signal-to-noise ratio (SNR) was calculated in the [0.0, 20.0] d"' frequency range. 

- The V and V-I magnitude averages (i.e., the zero frequency constant in their Fourier decompositions) are given in those cases in which the 
flux-to-magnitude conversion was reliable (see Sect.|2]and Fig.[8}. 

- Stars denote new detections, RRab, RRc and RRabB/RRcB stand for fundamental, first overtone and Blazhko variables. 

- Fourier decompositions of the RR Lyrae stars with magnitude-transformed light curves in V color are available electronically online at CDS. 
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Fig. 5. Standard V light curves of the 24 RR Lyrae stars for which the flux-to-magnitude conversion was performed. CVSGC 
numbers and periods in days are given in each case (see also Table [T] for additional data). All light curves were reconstructed by 
TFA. 



able number of stars, the discrepant position in the PLC relation whether the same photometric error (and not, e.g., non-cluster 

was coupled with an obviously discrepant amplitude of the light membership) can account for the outlier status of all the omitted 

curve. This is evidently introduced by the incorrect conversion stars. The correlation of the formal errors of {V - /)rf with the 

of their fluxes into magnitudes (see Sect 12. 11 ) due to the inaccu- (vertical) distance from the best-fitting PLC relation is exhibited 

racy of their RF photometry. However, it is important to justify in Fig. [8] Indeed, all omitted stars have also larger photometric 
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Fig. 6. Light curves expressed in arbitrary differential flux units of the 30 RR Lyrae stars that were not transformed into magnitude 
values due to the large zero-point errors associated with crowding. All plots are shown on the same ordinate scale. CVSGC numbers 
and periods in days are given in each case (see also Table[T]for additional data). All light curves were reconstructed by TEA. 



formal errors. Therefore, we think that the data clipping used in 
cleaning up the PLC relation is justified. On the other hand, the 
"good quality stars" (defined as such with a formal (V - /)rf er- 
ror of < 0.03 mag, located below the gray line in Fig. |8]l show 
various levels of agreement with the ridge line of the PLC, rang- 
ing from perfect match to near 3cr difference. This indicates that 
formal image errors can serve as a guide only when strong dis- 
crepancies are observed in applications involving photometric 
data derived from the lower quality images. The final standard 
deviation of the fit is 0.063 mag, slightly worse than the one de- 
rived from several clusters by KWOl. 

For the zero-point calibration of this PLC rela t ion we relied 
on the Baade-Wesselink (BW) analysis of lKovaci ( l2003b . From 
the 21 BW stars with V - I color indices we omitted AV Peg, 



because of its strongly discrepant position on the PLC relation. 
Six additional stars (X Ari, SW Dra, SS Leo, V445 Oph, BB Pup 
and W Tuc) could also be omitted but this would not lead to 
a significant change in the zero point of the BW PLC relation 
(-1.212 + 0.037 vs. -1.182 + 0.026 with the six stars omitted). 

By fixing the slope of the log Pq term as given by KWOl for 
W{VJ) = y-2.45(y-/), for the 24RRab andRRc stars of M53 
we get 



WiVJ) = -2.513 logPo + (15.100 + 0.013), 



(1) 



where the error corresponds to the Icr formal statistical error 
Similarly, for the 20 BW stars we get 



W(VJ)BVi = -2.513 logPo- (1-212 + 0.037). 



(2) 
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Fig. 8. Distance from the empirical PLC relation versus the er- 
ror of the color index computed from the formal errors of the ob- 
ject images on the reference frames. Gray horizontal line shows 
the near 3cr limit of the scatter around the regression line of the 
empirical relation of KWOl. Dashed vertical line indicates the 
cutoff error for the color index as explained in the text. 



Simple subtraction of the zero points of the above two equa- 
tions yields 16.312 + 0.039 for the dereddened distance modu- 
lus of the cluster. We note that this value is compatible with an 
LMC distance mo dulus of 18.5 mag ()Kovacsll2003l, see also the 
recent summary of Idi Benedettoll2008l) . Additional notes on the 
distance of the cluster will be given in Sect. l3.3l in relation to the 
metal content and various Mv-[Fe/H] relations. 



3.3. [Fe/H] distribution 

The application of the empirical for mula for the computation 
of the iron abundance of RRab stars (lJurcsik & Kovacslll996i 
hereafter JK96) has been proven to be quite successful in the 
past (Jurcsik & Kovacs 1996 [comparison for several GCs]; 
[Szekelv et alj 12007 [NGC362]). However, there is a growing 
concern o n the accuracy of this formu la at the low-metalhcity 
end (e.g., iNemed l2004l FNGC 50531: iGarcia Lugo et all l2007l 
[M15]; see however lKinemuchi et al ] |2008 l [Draco RRab vs. gi- 



ant metallicities]). Since M53 is also a low-metallicity cluster, 
and the available light curves are of high-quality for computing 
precise Fourier decompositions, here we can examine this prob- 
lem in more detail. 

In tracing back the meta llicity estimates der ived for this clus- 
ter, we found the works o f .Zinn & West! (11984 hereafter ZW84) 
and ISuntzeff et all d 19881 hereafter S88). The metallicities pre- 
sented in these papers are both based on calibrated spectral in- 
dices (similar to the AS index of Preston 1959 on which the 
JK96 formula is calibrated). The two studies, based on the anal- 
yses of several giant stars, yield very similar metallicities for 
M53 (-2.04 + 0.08 by ZW84 based on the Qjg parameter, and 
-2.09 + 0.09 and -2.20 + 0.33 by S88 based on the othk and 
mMg indices, respectively). High-disp ersion detailed abundanc e 
analysis on a single giant was made bv lPilachowski et akl (Il983[) . 
They obtained a value of - 1 .9. Additional early [Fe/H] estimates 




log(g) 

Fig. 9. Estimation of the gravity of t he M53 giants from col - 
ors by using the evolution models of IPietrinferni et all ( |2004|) . 
Overall iron abundances are shown at the left of the isochrones 
spreading from 10 to 13.5 Gyr (please note that the age depen- 
dence is highly degenerate). Horizontal bars denote the positions 
of the giants as listed in Table I A. 1] 



are also listed in their paper. Most of the values indicate a metal- 
licity lower than -1.8. 

To further confirm the low metallici ty of the giant s of M53, 
we employed the photometric method of lKovacsl (1200 8*) (see also 
Dekany et al. 2008). The method is based on the well-known ob- 
servation that the B - V color index depends fairly sensitively 
on the gravity and more importantly on the metallicity. Because 
the most relevant parameters determining the theoretical values 
of the broad-band colors are the effective temperature, gravity 
and overall metal content ([M/H], identical to [Fe/H], assuming 
solar-type heavy element distribution), we can estimate [Fe/H] 
by using BVI magnitudes, assuming that we can give a reason- 
able estimate on \ogg. For Cepheids and RR Lyrae stars this 
estimate is very accurate, since the period is in tight correla- 
tion with the gravity. For giants we can get estimates on the 
possible ranges of the gravity by employing evolutionary mod- 
els. We took the isochrones from the evolutionary models of 
Pietrinfer ni et al. ] (2004) with metallicities bracketing the esti- 
mated metallicity of -2.0 dex. Then, we checked the position 
of the selected giants (see Appendix |A]i with available three- 
color photometry on the logg —> (B - V), (V - I) diagrams. We 
then took the corresponding log g ranges and computed individ- 
ual me tallicities from the stellar atmosphere models of ICastellil 
(Il999h . 

The relation between log g and the color indices at various 
metallicity values and cluster ages is shown in Fig. |9] It is very 
fortunate that there is basically no age-dependence on this di- 
agram, so one can read off the corresponding logg ranges af- 
fected only by the adopted metallicity values. It turns out that 
in the color index - logg space of the giants studied, there is 
a very minute dependence of the derived metallicity on the as- 
sumed values of log g. (We note however, that for log g values 
exceeding the high-metallicity bounds by ~ 0.3 - 0.5 dex we 
could get a significant increase of up to 0.5 dex in [Fe/H].) With 
the log g ranges confined in this way, we computed more accu- 
rate individual photometric abundances by matching both color 
indices simultaneously. The derived values, together with other 
important data on the stars used, are given in Table lA.ll 
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For comparison with the values obtained on RRab stars by 
using the empirical formula of JK96, we plotted the metallicities 
for the two groups of stars in Fig. [TO] Error bars of the individ- 
ual RGB stars denote standard deviations, obtained by adding 
uncorrected Gaussian random noise with cr - 0.02 mag to their 
B—V and V- / color indices (the very small formal errors for the 
RRab stars are not shown). In spite of the large scatter, it is clear 
that the RGB stars have significantly lower overall metallicities. 
Indeed, the average metallicities with their formal 1 cr errors are 
-2.12 + 0.05 and -1.58 + 0.03, i.e., the difference is well above 
the 3cr level. We note that our average photometric metallicity 
shows fine agreement with the ones obtained by spectroscopic 
measurements, so the suspected systematic difference between 
the HB and RGB populations seems to be confirmed. However, 
we do not know at this point what is the source of this difference. 
Although we cannot exclude that there are substantial chemical 
inhomogeneities in some clusters (e.g., NGC 1851, see lLee et alj 
12009,) . we think that it is more probable that the JK96 formula 
is biased upward at lower metallicities due to the low num- 
ber of low-metallicity stars that entered in the calibration sam- 
ple (see a lso Schwarzenberg-Czernv & Kaluznv 1998; Kovacs 
l2002t iNemed 1200 4). To check if this is indeed the case, we 
would need either precise BVI photometry, or, even better, high- 
dispersion spectroscopy for the RR Lyrae stars. Unfortunately, 
there are no such data available. We made a humble attempt to 
utilize sor ne of the spora dic mu lticolor data for, e.g., M15, M68, 
and M 92 dCorwin et al J 120081: IWalkeiill994 iKovacs & JurcsikI 
1 19971 and references therein) but either the number of objects 
was low or the average magnitudes were too poor to derive con- 
sistent metallicities. 

Since there are ample amount of RRc star s in M53, we de- 
cided to employ the formula calibrated by Mor gan et al.l ( |2007b 
for an independent estimation of the metallicity of the RR Lyrae 
population. By us ing thei r Eg. (3) and transforming the val- 
ues to the scale of |Jurcsil3 d 19951) . we get an average [Fe/H] of 
-1.68 + 0.08 (error of the mean) for the sample of 15 stars (for 
two stars, data were taken from Kopacki 2000). Although this 
value confirms the one derived from the RRab population, we 
note the following: (i) the structure of the RRc light curves con- 
tain less distinct features than those of the RRab stars and even 
these are often washed out by noise: (ii) the calibration for the 
RRc stars is based on cluster variables that, because of the ab- 
sence of individual metallicity estimates, are assumed to have 
same metallicity as that of the host cluster (i.e., there are 12 in- 
dependent metallicity values used in the regression for the 106 
variables in the calibrating sample): because of the use of cluster 
variables, the full span of metallicity range is only [-1.0, -2.2], 
that is nearly a factor of two smaller than the one used by JK96 
for the Galactic field RRab stars. We presume that the good over- 
all agreement between the metallicities is primarily attributed to 
the strong dependence of [Fe/H] on the period for both types of 
stars. 

We made yet another test on the possible dependence of the 
derived metallicity on the chemical composition (namely on a- 
element enhancement, since the cluster is of low-metallicity). It 
turned out that the a-enhanced models shift the lower boundary 
of the log g range downward by some ~ 0.2 dex and resulted in 
a small (i.e., < 0.1 dex) and apparently non-systematic change 
in the metallicity [M/H]. We recall that if the solution of the dis- 
crepancy between the JK96 and RGB metallicities lays in the 
non-solar chemical composition, then we would expect a dra- 
matic increase in the overall RGB metallicity [M/H] estimated 
by the above method from the a-enhanced models, since then 
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object 

Fig. 10. Photometrically derived overall iron abundances of 12 
RGB stars of M53 (black diamonds) versus iron abundances of 
16 RRab stars computed from the V light curves by the JK96 
form ula (gray circle s [data from this paper] and squares [data 
from Kopacki '2000*]). Horizontal lines show the correspond- 
ing averages (for comparison, we also show the value given by 
IZinn & Westill984i bv dashed line). Formal Icr errors are given 
for the means. For star IDs we refer to Tables [Tland lA. II 

the relative abundance of iron (which the JK96 formula is cali- 
brated to) would be lower in these models. 

Following the recommendation of the referee, based on the 
above results we can check how the cluster's distance, as derived 
from the PLC relation, relates to the one estimated from the more 
tradi tional (but less ac curate, due to evolutionary effects, see, 
e.g.. ICassisi et alj 12004) Mv-[Fe/H] relation calibrated by the 
same BW distance scale. Relying on the same BW sample as 
in Sect. 13.21 and using the same spectroscopic metallicities as 
in ' Dekanv et alj (|2008|) we get the following linear calibration 
between the absolute V magnitude and the iron abundance: 

Mv = (0.24 ± 0.06) [Fe/H]jK96 + (0.88 + 0.06) . (3) 

As expected, the standard deviation of the above fit is huge 
(0.17 mag as compared to 0.063 mag of the PLC relation - see 
Sect. 13.21 ). Although many variants of Mv-[Fe/H] calibrations 
can be found throughout the literature, the slope of the above for- 
mula turns out to be very close to the value of 0.23 ±0 04, fa vored 
by recent reviews in this subject (see, e.g.. ICatelanI 12009. and 
references therein). The mean V magnitude of M53 RR Lyrae 
stars is equal to 16.85 + 0.01 mag, based on the magnitude- 
transformed light curves of 24 stars (see Table[Tli. Using the two 
different metallicity values derived for the RR Lyrae and red gi- 
ant stai-s, from Eq. (O we get 16.35 + 0.11 and 16.48 + 0.14, 
respectively, for the distance modulus of the cluster. Both are in 
agreement with the PLC result in Sect. [32] within the errors, but 
with an apparent preference for higher metallicity. 

4. Short-period variables 

We Fourier-analyzed the photometric time-series of all identified 
objects in the [0, 50] d ' frequency range to search for short peri- 
odic variables, in particular for SX Phoenicis and high-amplitude 
Delta Scuti stars. With the aid of TEA, we considerably increased 
the reliability of detection of short-period, low-amplitude signals 
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Fig. 11. Left panels: TFA-reconstructed and folded light curves 
of the short periodic stars. Right panels: DFT amplitude spectra 
of the corresponding time-series. The separation between ordi- 
nate ticks denotes the same flux difference in each plot. IDs and 
periods are given in each case. See Table |3] for additional prop- 
erties of these variables. 
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Fig. 12. TFA-reconstructed light curves and DFT amplitude 
spectra of V77 and V8 1 . Labeling is the same as in Fig. [TT] 



by filtering out low-frequency trends from the data whose am- 
plitudes would have overwhelmed the intrinsic light variations 
of these objects. 

All of the short-period stars previously known in M53 were 
discovered by lJeon et al ] (l2003h who used a 1.8-m telescope and 
PSF photometry. All of these objects are of SX Phe-type lo- 
cated in the blue straggler region of the color-magnitude dia- 
gram (CMD). We successfully identified SXPl, SXP2, SXP6, 
and SXP7 among these eight objects in our data (they are identi- 
cal to our V7 3, V74, V75, and V 76, respectively). From the data 
published by Jeon et al.1 (|2003|) . SXP4, SXP5, and SXP8 have 
too low amplitudes to be identified in our data. The situation is 
different for SXP3 w hich has a reasonably large amplitude of 
0.1 mag according to iJeon et al.l (2003) but it is rather close to 
V66 and V67 and the limited resolution of our CCD chip did not 
enable us to separate their light variations. 

In addition to the previously known SX Phe stars, we dis- 
covered seven new objects (signed as V77 through V83) that 
have significant high-frequency signals. Except for two of them, 
these faint objects were heavily crowded or merged with brighter 
stars, therefore, we were not able to obtain accurate RF mag- 
nitudes for them to perform the flux-to-magnitude conversion. 
To investigate if these objects were located in the blue strag- 
gler star (BSS) region, i.e. they could be securely classified as 
SX Phe stars, we sought for supplementary information in the 
literature. We cross-correlated the positions of the 7 new objects 
with th e single-epoch data-base of M53 BSSs bv iBeccari et al.l 
(2008a, via private comm.) and found a match with sub-pixel 
accuracy in three cases, namely V77, V78, and V79, which sug- 
gests that these variables are most probably of SX Phe type. 
For V77 {V = 18.478, V - / = 0.433) and V81 {V = 17.454, 
V -I - 0.815), we were able to perform the conversion of fluxes 
into magnitudes. The average magnitudes of V77 verifies its lo- 
cation in the BSS region, while V81 is l ocated exactly on th e 
lower red giant branch (see the CMD of iBeccari et al.ll2008al) . 
We are not sure whether there is a blending eff'ect of two stars 
in this case or V81 (which we identified as USNO-Bl.O 1081- 
0245849) is indeed a peculiar variable. We note however that 
there is a suspected BSS (USNO-Bl.O 1081-0245846) near V81 
which might account for this variability, although we did not de- 
tect the corresponding signal in its time-series. 

The type of the remaining three new short-period variables 
(V80, V82, and V83) remains unknown, but their blue strag- 
gler nature cannot be ruled out (for example, we note that the 
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Table 3. Properties of the short period stars 



ID 


a [hms] 


S [dms] 


dl"] 


/[d-1] 


SNR 


V73 


13:13:03.1 


18 


09:26.1 


119.9 


14.26581 


10.2 


V74 


13:12:49.7 


18 


07:26.3 


181.9 


22.03949 


13.7 


V75 


13:13:09.4 


18 


09:39.5 


203.9 


22.60022 


11.8 


V76 


13:13:05.0 


18 


08:35.3 


168.4 


24.10224 


10.4 


V77* 


13:13:20.8 


18 


15:34.9 


488.3 


13.00995 


23.0 


V78* 


13:12:49.8 


18 


08:56.0 


106.9 


22.25488 


12.7 


V79* 


13:12:46.6 


18 


11:36.9 


150.6 


21.59184 


9.7 


V80* 


13:12:57.4 


18 


10:13.8 


31.3 


14.83029 


14.2 


V81* 


13:13:02.7 


18 


06:29.4 


244.5 


14.01180 


9.5 


V82* 


13:12:56.5 


18 


13:10.0 


181.1 


45.33103 


8.4 


V83* 


13:12:50.1 


18 


07:42.9 


164.5 


8.01907 


8.2 



Notes: New di.scoveries are marked with asteri.sks. The .signal-to-noi.se ratio 
(SNR) was calculated in the [0, 50] d"' frequency range. See Table[T]for other 
notations. 



previously known va riable V73 was not identified as a BSS by 
iBeccari et al.ll2008 a). The properties of all 1 1 short-period vari- 
ables detected in this study are given in Table[3] The light curves 
and amplitude spectra of the objects are shown in Figures [TTHT2l 



5. Long-period variables 

Most of the stars on the upper part of the giant branch of M53 
show some kind of low-amplitude light variations on time-scales 
of greater than 10 d. We identified V65, V66, V67, V69, and 
V70 (see the CVSGC for their positions) among the previously 
known variable red giants, all of them showing a large power 
excess at very low frequencies (~ 0.01 d"'), implying long-term 
irregular variations. 

At the same time, we found three new long-period variables 
(LPVs) among the red giant stars (signed as V84, V85, and V86). 
The properties of these new objects are summarized in Table |4] 
and their light curves are shown in Fig. [T3] Variable V86 was 
saturated in the / band, therefore we give only an approximate V 
magnitude in Table|4]as a reference. V84 and V86 have very sim- 
ilar periods. We excluded the possibility of detecting the same 
systematics in these objects by using a large variety of template 
sets (see Sect. 12.2b and also by verifying that the two light vari- 
ations are not in the same phase. 

Recent investigations based mainly on OGLE-II data have 
showed that most of the pulsating red giants below the tip of the 
red giant branch ( TRGB) are among the first ascent RGB stars 
(see e.g.. lKiss & B edding 2003). The majority of these stars are 
sitting on parallel ridges in the period-luminosity (PL) plane. We 
inspected whether the above two LPVs with V - I color indices 
do match with a P L ridge on the log P-W esenheit index plane 
of LMC giants by ISoszvnski et all ( l2005b . assuming an LMC 
distance modulus of 18.5. We found that V84 is below the faint 
limit of the LMC ridges, while V85 shows a fine match with 
rid ge A/R3 wh ich was associated with third-overtone pulsators 



ge A/K3 wn 
bv lWoodldSoOOl) 



6. Conclusions 

Globular cluster M53 was put on our observing plan with the 
prime goal of extending the available data-base on its variable 
stars. This is the first time-series photometric study covering the 
whole cluster up to ~ 14' (for a comparison, the tidal radius of 
the cluster is ~ 22', according to iHarrisll 19961) . We analyzed the 
whole set of V, I photometric time-series data acquired on 3048 



objects in two seasons in 2007-2008. Our goal was to establish 
a reasonably deep variable star inventory by considering all pos- 
sible variable types that may occur throughout the Hertzsprung- 
Russell diagram. Because of the limited time span of the obser- 
vations, the search was limited to periods < 100 days. Similarly, 
the lower bound of the brightness level was limited by the tele- 
scope aperture and overall sky conditions to V < 21 mag, in- 
cluding the blue straggler region a nd ending ~ 1 m ag below the 
main sequence turn-off point (e.g.. lRev et al.ll998l) . In analyzing 
the data we used the OIS method for image processing and TEA 
for time-series post-processing. 

Altogether 12 new variables were identified, most of 
them are of short periodic ones (three are of SX Phe-type). 
Interestingly, and somewhat surprisingly, except for the three 
Blazhko variables, we did not find any multiperiodic variables, 
although they could be quite abundant among RR Lyrae and 
SX Phe stars (iKovacsl 120011: iGilhland et all Il998h . We were 
also unsuccessful in finding eclipsing binaries, although some 
10% of the blue stragglers are thought to form binary sys- 
tems (Beccar i et al.ll2008bl) . We note that in NGC 5466 , a sim- 
ilarly low metallicity cluster. [Arellano Ferro et aL I (120081) identi- 
fied th ree eclipsing variables discovered earlier bv iMateo et alj 
(11990 ). They are all in the BSS region. It may well be that the rel- 
atively high scatter of our data at their low brightness level dis- 
abled us to detect these types of binaries. It is interesting to note 
that we did not find variables slightly more luminous than the 
horizontal branch near the red giants. In a simil ar analysis of a 



sample of the Large Magellanic Cloud field stars, Szulag vi et al 
(2009) found several such stars (see also Arellano Ferr o et al 



2008 for two such stars, one is right on the red giant branch of 
NGC 5466). We note however, that we have four short periodic 
stars with periods reminiscent of the SX Phe stars. Although, ex- 
cept for one object, we do not have color information on these, 
based on thek very low periods of less than 0.13 days, we sus- 
pect that all of them are actually SX Phe variables. 

In analyzing the variables identified in the present data, we 
focused on the RR Lyrae stars. Two important parameters (av- 
erage distance and metallicity) that are relevant for the physi- 
cal properties of the cluster were derived from the light curves 
and average colors. We have shown that the period-luminosity- 
color (PLC) relation spanned by the RRab stars as observed 
in other globular clusters, is also very well exhibited in M53. 
Furthermore, it was also demonstrated for the first time on em- 
pirical basis that RRc (first overtone) RR Lyrae stars, after proper 
period shift, follow the same linear log Po-W(V - I) relation 
as those of the RRab (fundamental mode) s tars. By using the 
Baade-Wesselink calibration of this relation (' Kovacsir2003 ^. we 
derived the distance modulus of 16.3 1 + 0.04 mag for the cluster. 

By using the iron ab undance-Fourier parameter relation of 
lJurcsik&Kovacsl ( 119961) . we derived an overall iron abundance 
of -1.58 + 0.03 dex for the RRab stars. This is significantly 
higher than the value of -2.1 computed from low-dispersion 
spectra for the red giant stars. The lower metallicity of the gi- 
ants has been confirmed by our estimate based on BVI pho- 
tometry, stellar atmosphere and evolution models. Such a dis- 
crepancy of 0.3-0.5 dex, seems to be present also in other low- 
metallicity clusters. For example, for the RR ab stars of M68 we 
find [Fe/H]= -1.8 dJurcsik & Kovacsll 19961). whereas the high- 
dispersion spectroscopic data on giants by Lee et al.l (|2005[) sug- 
gest a value of -2.2 dex (both values are on the scale of 'Jurcsi^ 
1995). Similarly, for the RRab stars of NGC 5053 Nemec (20041 
finds [Fe /H]= -1 .7, wh ereas the low-dispersion spectral indices 
of iSuntzeff et al.l (Il988b yield -2.3 dex for the giants. 
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ID 


a [hms] 


6 [dms] 


d["] 


Pld] 


SNR 


V 


V-1 


V84 


13:12:36.2 


18:07:32.3 


306.6 


22.4 


15.5 


14.766 


1.175 


V85 


13:12:50.7 


18:10:39.7 


69.7 


19.8 


19.2 


13.962 


1.357 


V86 


13:12:52.7 


18:10:28.0 


39.4 


22.2 


12.9 


~ 14.2 





Notes: For variable V75, only instrumental V magnitude average is shown because the star was saturated in the / band. The signal-to-noise ratio (SNR) was 
calculated in the [0, 20] d"' frequency range. 

V86 ' P = 22.2cl" 

0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

phase phase phase 

Fig. 13. V light curves of the 3 new long-period red giants folded with their periods. Star IDs, periods, and sine curves fitted to the 
data are also shown. 




Our current understanding of this discrepancy lays in the 
lack of enough number of low-metallicity RRab stars in the orig- 
inal calibration set of Jurcsik & Kovacs (1996). Unfort unately, 
from the works of ISuntzeff et al.l (1 1994)) and that of iLavdenl 
( Il994i) there were no new efforts for obtaining metal abun- 
dances for large number of field stars. A reliable re-calibration 
of the formula cannot be made without a new and homoge- 
neous sample. As an example for the need of new spectro- 
scopi c survey, we mention the case of UU Boo dJurcsik et al.l 
120081) . for which star the current Fourier formula yields [Fe/H] - 
- 1 . 17, whereas the spec troscopic observations of lLavdenI d 19941) 
and iKinman & Carrettal (1992) resulted in values of -1.64 and 
-1 .00, respectively. Although other tests based on new and accu- 
rate measurements of non-Blazhko stars sh ow remarkable ag ree- 
ment with earlier spectroscopic values (e.g.. lKun et al.l2008 ). the 
above example strongly supports the need of revisiting Galactic 
field RR Lyrae stars by high-dispersion spectroscopic surveys. 
Within the same framework, a more thorough discussion of the 
RRab - giant metallicity discrepancy would require a systematic 
high-dispersion spectroscopic study of the RR Lyrae and other 
non-variable populations of globular clusters. Using the same 
type of measurements of all cluster members selected and an- 
alyzing the chemical composition in the same way would lead 
to a uniform treatment of their composition and, in particular, 
to the avoidance of the continuous confusion about the differ- 
ent zero points and scales (i.e., Zinn & West 1984 vs. high- 
dispersion values). Importance of the more sophisticated treat- 
ment of the metallicity issue is underlined by current works 
on detailed a bundance analyses of globular clust er stars (e.g.. 



Lee et al. 2009) and on individual field stars (e.g.. lPreston et al.] 
2006.; Wallerstein et al...2009,) . 
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Appendix A: Red giant stars 

For the estimation of [Fe/H] using the photometric method de- 
scribed in Sect. 13.31 we needed accurate B - V and V - I color 
indices for a number of red giant stars in M53. Since our obser- 
vations covered on ly the V and / b ands, we took the B-V values 
from the catalog of Re v et"an ( ll99 8). We selected a large number 
of stars from this catalog with V and B-V values consistent with 
an RGB evolution ary state accord ing to the color-magnitude dia- 
gram published bv lRev et alj (Il998i) . We mapped the instrumen- 
tal pixel positions of the catalog onto our reference system by 
using the positions of RR Lyrae stars common in Rey's catalog 
and our data-base. This transformation yielded positions with an 
accuracy of a few arcseconds which was sufficient to unambigu- 
ously identify these bright stars on our frames by eye. As a fur- 
ther selection procedure we rejected the saturated stars and kept 
only those which were well outside the core and therefore had 
reasonably low photometric errors on the RFs. After the flux-to- 
magnitude transformation, we computed the simple V - I mag- 
nitude averages for these objects. We note that almost all RGB 
objects in our sample showed some kind of low amplitude ir- 
regular light variations on long timescales. Therefore we further 
rejected those objects for which th e difference was too high be- 
tween our V magnitude and that of lRev etalJ(ll998h . In our final 
sample of 12 RGB stars the standard deviation in the differences 
of the two V magnitudes is 0.02 mag. Table lA. ll lists catalog IDs, 
coordinates, distances from the cluster center, magnitudes, and 
estimated [Fe/H] values for these 12 red giants. 
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ID 


a [hms] 


d [dms] 


d["] 


V 


B-V 


V-I 


[Fe/H]Bv/ 


20 


13:12:56.9 


18:08:41.9 


91.2 


14.04 


1.30 


1.35 


-1.91+0.25 


23 


13:13:03.7 


18:09:38.7 


123.9 


14.07 


1.32 


1.38 


-2.09 + 0.25 


27 


13:12:37.4 


18:08:23.4 


276.5 


14.08 


1.32 


1.38 


-2.10 ±0.25 


29 


13:12:50.4 


18:08:53.7 


102.6 


14.11 


1.29 


1.37 


-2.24 ± 0.22 


40 


13:12:48.6 


18:14:15.7 


263.5 


14.34 


1.20 


1.28 


-2.04 + 0.25 


51 


13:13:17.4 


18:14:46.0 


419.4 


14.56 


1.08 


1.21 


-2.35+0.20 


53 


13:12:48.2 


18:12:46.6 


186.3 


14.59 


1.08 


1.20 


-2.31 ±0.20 


54 


13:13:05.3 


18:14:50.8 


315.7 


14.62 


1.07 


1.20 


-2.32 ±0.21 


61 


13:12:55.6 


18:13:13.8 


184.2 


14.74 


1.05 


1.15 


-1.89 ±0.27 


67 


13:12:38.3 


18:11:04.7 


248.0 


14.80 


0.99 


1.12 


-1.98 ±0.27 


68 


13:13:02.5 


18:13:23.2 


219.2 


14.83 


1.03 


1.15 


-1.98 ±0.26 


72 


13:12:47.9 


18:06:32.2 


241.5 


14.88 


1.04 


1.14 


-1.85 ±0.26 



Note: ID numbers and B-V values were taken from the catalog of lRev et"an (Il998h . Photometric metallicities [Fe/HJay; have been derived as 
described in Sect. 13.31 



